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ABSTRACT 


A 

It  Is  found  that  , 


A f  1°' 


i  te  perturbation  of  the  ion  orbits 


leads  to  a  nonlinear  frequency  shift  that  reduces  the  mode  frequency 
and  has  a  weak  stabilizing  effect  on  the  lower-hybrid  drift  instability. 
This  result  is  obtained  from  a  self-consistent  solution  of  the  Vlasov- 
Poisson  equations  using  perturbation  theory  in  which  the  nonlinear  die¬ 
lectric  function  and  the  nonlinear  temporal  evolution  of  a  single  unstable 
mode  In  the  low  drift  velocity  regime  are  calculated  analytical ly. 
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I.  INTRODUCTION 

The  linear  theory  of  the  lower-hybrid  drift  instability  is  well 
understood  and  has  been  discussed  in  detail  by  Davidson  et  a  1 . 1  When 
the  amplitude  of  the  wave  is  small  but  finite  after  a  time  equal  to  many 
multiples  of  the  growth  time,  further  evolution  will  be  different 
from  exponential  growth  at  very  small  amplitudes.  In  order  to  analyze 
this,  the  nonlinear  dielectric  response  function  and  the  nonlinear  temporal 
evolution  of  a  single  unstable  mode  are  derived  self-consistently  by  using 
perturbation  theory  to  solve  the  Vlasov  and  Poisson  equations.  The  single¬ 
mode  approximation  is  valid  for  the  instability  when  the  plasma  parameters 
are  close  to  those  for  linear  marginal  stability.  This  requires  vE<<vt; 

(the  low  drift  velocity  regime)  for  the  lower  hybrid  drift  instability, 

2 

where  v  ■  c£/B  is  the  equilibrium  ExB  drift  velocity  and  v  . -T./M  is 
E  tii 

the  ion  thermal  speed.  M  and  T.  are  the  ion  mass  and  temperature,  res¬ 
pectively.  Similar  single  mode  studies  include:  modulation  of  the  Lang- 

2  1 

muir  wave  due  to  weak  nonlinearity  9  ;  nonlinear  evolution  of  drift- 

4-7  8 

cyclotron  and  drift-cone  instabilities  both  in  theory  '  and  simulation. 

In  this  paper  we  demonstrate  that  a  finite  amplitude  perturbation  of 
the  ion  orbits  develops  during  growth  which  leads  to  a  weakly  stabilizing  non¬ 
linear  shift  of  the  lower-hybrid  drift  mode  frequency.  The  largest  shifts  in 
the  mode  frequency  occur  for  modes  with  wavelengths  much  longer  than  that  of 
the  most  unstable  mode  (k  «k  ) ,  at  which  wavenumbers  the  lower-hybrid  drift 
instability  converts  Into  the  drift  cyclotron  Instability. 

For  simplicity,  we  use  a  one -d 1  mans  1  one  1  slab  configuration  shown  in 
Fig.  1;  wave  propagation  is  in  the  x  direction;  the  magnetic  field  is  uni¬ 
form  and  In  the  z  direction;  the  density  gradient  is  in  the  y  direction. 


r  it  -  dn  i  Wwarn  ■  "i  i  m  imfcl 
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The  ions  are  Created  as  unmagneti zed  because  Che  wave  frequency  and  growth 
rate  are  much  greater  than  the  ion  cyclotron  frequency.  The  ions  are  in 
force  balance;  that  is,  the  force  due  to  the  ion  pressure  gradient  in  y  can¬ 
cels  that  of  the  equilibrium  electric  field.  The  lower-hybrid-dri  ft  instab¬ 
ility  is  analyzed  in  the  electrostatic  limit;  electromagnetic  effects  are 
assumed  to  be  small. 

In  Sec.  II,  the  nonlinear  dielectric  response  is  calculated  by 
solving  the  coupled  Vlasov-Poi sson  equations,  and  a  nonlinear  dispersion  rela¬ 
tion  is  obtained.  Section  III  is  devoted  to  a  derivation  of  the  time  evolu¬ 
tion  of  the  lower-hybrid  drift  instability.  The  field  energy  level  at  which 
saturation  might  occur  and  the  frequency  shift  due  to  the  finite  amplitude 
of  the  wave  are  also  determined.  Finally,  conclusions  and  a  comparative  dis¬ 
cussion  of  several  saturation  mechanisms  are  given  in  Sec.  IV. 

II.  DERIVATION  OF  THE  NONUNEAR  DIELECTRIC  FUNCTION 

2 

We  follow  the  method  of  reductive  perturbation  theory.  Ue 
assume  that  the  distribution  functions  FS(y,v,t)  for  species  s  and  the 
electric  potential  <*(x,t)  can  be  expanded  as 

•  P*(y.v)  ♦  £  cnP* (y ,v,t)ein8  ♦  c.c.  (1) 

n*  1 

and 

#(x,t)  •  2  (x,t)e'n0  ♦  c.c.  (2) 

n-1 

where 

*!(y.v,t)  •  2  F* , (y . V . t )  ,  n  •  0, 1 , . . .  (3) 

"  j-0  nj 


(4) 
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♦„(x,t)  -  E  .\jtx.t)  .  '.2 - 

j-0  J 

0  =  kx  -  ut  •  (5) 

The  small  parameters  is  e,  which  is  on  the  order  of  ed/T.<x1.  k  and  u 
are  the  wave  number  and  frequency  of  a  single  mode.  Quasi  linear  analy¬ 
sis  indicates  that  current  relaxation  (the  relative  drift  between 

q 

the  electrons  and  ions  goes  to  zero)  can  cause  saturation  for  vg<vtp 
However,  the  effect  of  current  relaxation  is  small  for  v£4evc;  •  There¬ 
fore,  we  shall  specialize  to  the  case  v£<Xvc;  an<*  treat  the  density  gra¬ 
dient  and  Vg  as  constant  in  our  derivation.  The  distribution  function 
F^j (y,v,t)  can  then  be  expressed  as 

-  n0Cy)fnj(v,t)  .  (6) 

The  Poisson  equation  of  the  system  Is 

-  7%  -  4»noe  j dv  (f*  -  f*)  .  (7) 

Substituting  Eqs.  0}  and  (2)  into  Eq.  (7)  yields 

(nk)2in  -  4imoe  J  dv  ^  -  ?*)  .  (8) 

Since  the  characteristic  frequency  of  the  lower-hybrid  drift  instability 
Is  much  less  than  the  electron  plasma  and  cyclotron  frequencies,  it  is 
assumed  that  electrons  respond  to  the  wave  linearly,  i.e., 

■  *  x.("fc,n«) (nk)2i 

•  n 


kn n  e  f  dv  f* 

A  I  e  A 


(9) 
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where  xfi  is  the  linear  electron  susceptibility.  To  justify  further  the 
assumption  of  a  linear  electron  response,  we  also  assume  a  value  of 
zero  for  the  plasma  beta  (plasma  pressures  magneti c  pressure)  and  » 0  , 
so  that  electron  resonance  broadening  can  be  neglected.  For  finite  elec¬ 
tron  temperature  and  plasma  beta,  electron  resonance  broadening  can  stab- 

•  1  J  _  ,L  .11  .  *  ^  10,11  |  |  *  p  [  —  \  /  I  \  J  /  ft  \  p  /  Q  . 


i 1 ize  the  instabi 1 i ty. 


duces  to 


Using  Eqs .  (3),  (*♦)  and  (9),  Eqs.  (8)  re- 


(1  +  (nk.n*u)  ]  (nk)  .  «  kvn  e  f  f 1  .  dv 

e  nj  o  J  nj 


The  Vlasov  equation  for  the  ion  distribution  in  one  dimension  is 


dt  3x  M  3x  3v 


which  can  be  rewritten  as 


hL  =  /j_  +  v ±\f>  .  th.uL 

dt  \  3t  3x  /  H  3x  3v  ’ 


so  that  f  can  be  integrated  over  the  characteristics  corresponding  to 
the  unperturbed  orbits,  i.e.,  x»vt.  We  obtain,  as  usual, 


f1  -  f1 


00  « 


•  r  hisL^ 


r  h 

L  3* 


dt*  . 


Expansion  of  Eq.  (13)  in  the  series  given  by  Eqs.  (1)  and  (2)  yields 


'  1 
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Convergence  of  the  series  solution  for  the  perturbed  distribution  function 
requires  that  the  electric  potential  cause  only  a  small  perturbation  to 
the  unperturbed  orbits.  Hence,  ion  trapping  is  excluded  by  this  assumption. 

One  should  also  note  that  the  superscript  i  for  ions  has  been  dropped 
from  Eq.  (14).  Equating  coefficients  of  exp(in0)  at  the  same  order  in  e  in 
Eq.  (14)  yields 

fQ1  -  0  ,  (15) 

and 

f  .  !liO  .  8  fOO/3V  (16) 

10  M  v  -  V 

where  '/•w/k  and  f^g  is  the  equilibrium  ion  distribution  function  in  the 
absence  of  perturbation.  We  then  substitute  Eq.  (16)  into  Eq.  (10)  to  obtain 
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[*  ♦  Xe(k#tt)] 


k% 


4irn  e^ 

o 


10 


10 


f 


>W»» 

V  “  V 


dv 


k  X,  (k,u»)«10  , 


07) 


which  yields  the  linear  dielectric  function  D(k,iu), 


0 (k , tu)  -  1  +  xe(k,u)  +  x,  (k,u)  . 


08) 


We  now  proceed  to  calculate  the  second  order  components.  Assuming 
oj»ur+i«  and  <S-*0,  and  equating  coefficients  of  the  constant  terms  in  Eq. 
04)  gives 


09) 


Substituting^  by  using  Eq.  (16),  we  obtain  the  quasi  linear  modification 
to  the  distribution  function, 


The  component  f^  results  f row  the  second  order  terms  of  Eq.  (14)  for  exp(id), 


f 


11 


e*ii  3foo/,3v 

M  v-  V 


(21) 


r.. 


Then,  following  the  steps  used  in  deriving  Eq.  (17)  and  Eq.  (18),  we 
obtain  from  Eq.  ( 1 0) 


D(k,u»)<j>. 


The  components  f..  and 


are  uncoupled  from  the  lower  order  compo¬ 


nents  and  are  irrelevant  to  the  rest  of  the  calculation.  Therefore,  we 
set  d>  1 1  *  f  1 ,  *  0  in  the  rest  of  the  analysis. 

Equating  the  coeff i cients  of  the  second  harmonic  terms  (2ji,2k) 
in  Eq.  (14)  give 


i/*.oY  ■ 


!VH\  3foo/3v 


M  /  v-V  3v 


The  first  term  appearing  on  the  right  side  of  Eq.  (23)  is  the  modification 
due  to  the  unshielded  second  harmonic  oscillation  of  a  single  wave,  and 
the  second  term  represents  the  shielded  effect.  Similarly,  Eqs.  (10)  and 
(23)  yield 


3k  0(2k,2u) 


'*10  f  _ 3_ 

M  ^  v-V  3v 


3  f0()/3v 


For  the  third  order  component  f^,  Eq.  (14)  gives 


v-V 


3foo/3v 


"V3* 
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Using  f1Q  ,  fQ2  and  f^Q  from  Eqs .  (16),  (20)  and  (23),  we  find  that 


Substituting  Eq.  (26)  into  Eq.  (10)  for  $12,  combining  to  obtain  ^  ^  »<f>  1  1 2 , 

and  then  using  Eqs .  (2),  (17)  and  (22)  yields  the  nonlinear  dispersion  relation 


+®<e4)  • 


(27) 


1 1 


The  first  term  on  the  right  side  is  the  nonlinear  coupling  of  the 
potential  with  the  quasi  linear  perturbation.  The  last  two  terms  arise  from 
the  nonlinear  coupling  of  second  harmonic  variations  in  f'  and  $  with 
the  perturbations  at  the  fundamental,  with  and  without  plasma  shielding 
effects . 

We  define  W(z)  by 


W(z) 


2  r3W3v 


dv 


(28) 


where  *u>/kvt.  .  The  quasi  linear  term  gives 


/  -1 


a2  /3foo/3v 


v-V  3v2\(v-V)2 


_J _  d  w(z) 

12v‘i  to3 


The  unshielded  second  harmonic  effect  becomes 


and  the  term  associated  with  the  shielded  second  harmonic  oscillation  is 


u  - 


8y  using  Eqs.  (29),  (30)  and  (31),  we  rewrite  Eq.  (27)  as 

*|2¥  (32) 

where  i|»«e$j/T.  and  X ^  is  the  ion  Oebye  length.  The  quasi  linear  term 
and  the  unshielded  second  harmonic  effect  are  combined  in  the  first  term 
on  the  right  side  of  Eq.  (32). 

III.  NONLINEAR  EVOLUTION  AND  FREQUENCY  SHIFT 

In  this  section,  we  estimate  the  field  energy  at  saturation 
caused  by  a  finite  nonlinear  frequency  shift  by  solving  the  nonlinear 
dispersion  relation.  When  the  wave  amplitude  is  very  small,  Eq.  (32) 
reduces  to  the  usual  linear  dispersion  relation 

0(k,u»)  -  0R(k,a.)  +  iOj  (k,w)  .  (33) 

Let  us  examine  Eq.  (33)  in  the  lower  drift  velocity  regime  characterized  by 
Ir/w r [  «  I  ,  V  and  v£  «  v^.  ,  (34) 

where  and  y  are  the  real  and  imaginary  parts  of  frequency,  respect¬ 

ively.  The  dielectric  function  for  cold  electrons  and  Maxwellian  ions  is 


D(k,wH 


1 


I6k2xj 


d\<(z) 


1 


2k2X20(2k,2w) 


expressed  as 


When  the  amplitude  of  the  wave  is  small  but  finite,  we  expand  Eq. 


(32)  around  <u  by  replacing  u»  with  u  +  i  ( 3/3 1 )  and  use  Eqs.  (35) 
o  o 

and  (36) ,  to  obtain 


i[0  +  ia)  ~  -y]  ^  *  (A  +i  B)  j  | 2  ^ 


(41) 


where 


a  =  (30  j/3u)^'(30^/3(i))o  %  -y/u»o  ,  (42) 

A  .  (olo'kVE >2  A  ,  I 

8kv£  \3  AfoUk.J^) 


and 


B 


6 


By  using  Eqs.  (36)  and  (33) »  we  get 


(44) 


(45) 


Substituting  Eq.  (45)  into  Eqs.  (43)  and  (44),  and  using  Eqs.  (29)  through 
(32),  the  relative  strengths  of  the  nonlinear  contributions  from  the  quasi - 
linear  modification  (ql),  the  bare  second  harmonic  oscillation  (b)  and 
its  shielded  effect  (s)  are  given  as  Che  ratios, 


-  15  - 


'  :  A<2k,2“0)» 


;«> 


3iVJ 


ql 


(2k. 2*  ). 

O  0 


(2k.2u  ) 

o  s 


20  :  -15 


(1*7) 


For  the  most  unstable  mode,  k«k  ,  the  three  nonlinear  contributions  to 

nt 

the  ion  distribution  function  are  comparable.  For  k/k  «1  the  shielded 

fll 

second  harmonic  contribution  dominates  the  nonlinear  modification  of  the 

ion  response.  This  is  because  the  second  harmonic  perturbation  is  very 

close  to  satisfying  the  linear  dispersion  relation  for  k«  k  ,  i.e.,  D(2k,2w  ) 

mo 

nearly  vanishes. 

In  order  to  obtain  the  time  evolution  of  the  wave  amplitude  and 
frequency  shift,  we  define  t It  •  r  exp(-is),  where  both  r  and  s  are  real.  Eq. 
(42)  becomes 


ar  -  rs 


(48) 


and 

•  •  a  3 

r  +  ars  -  yr  ■  Br 

Eliminating  r,  we  obtain 

•  ay  A  ♦  aB  2 

S  -  - 4-y  «•  - J  r 

I  ♦a  1+o 

where  the  first  term  is  Che  linear  correction  to  the  frequency  in  the 
presence  of  growth,  and  the  second  term  is  the  nonlinear  frequency  shift 
which  grows  in  time  with  r  (i.e.,  |ed^/T.j  ).  Eliminating  rs ,  we  obtain 


(49) 


(50) 


Y  B  -  oA 
1  +a2  1+a2 


Integration  of  Eq.  (51)  yields 


r 


2 


2  ce2Yt/^  +a2^ 

r’  ~ 


where 


(5D 


(52) 


(53) 


is  the  field  energy  level  at  saturation. 

From  Eqs.  (42)  through  (47),  it  is  obvious  that  the  quasi  linear 

effect  and  the  shielded  second  harmonic  oscillation  stabilize  a  single 

lower-hybrid  wave,  i.e.,  they  nonlinearly  reduce  growth.  The  nonlinearity 

due  to  the  bare  second  harmonic  oscillation  enhances  growth  and  raises  the 

saturation  level.  If  r  »  r(t*0)  •  r  ,  Eq.  (52)  gives 

•  o 


(54) 


and  Eq.  (50)  becomes 


T 


Jll 

,2  1  ♦  a 
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We  identify  the  second  term  on  the  right  side  of  Eq.  (55)  as  a  nonlinear 
shift  of  the  mode  frequency. 

IV.  CONCLUSION 

A  nonlinear  dispersion  relation  for  the  lower-hybrid  drift  instab¬ 
ility  was  derived.  We  obtained  the  saturation  field  energy  and  the  nonlinea 
frequency  shift  at  saturation.  With  use  of  Eqs.  (35)  through  (44),  Eq.  (53) 
gives 


at  saturation  in  the  limit  Vg/vtf«1. 

In  Fig.  2  we  present  !54t»  the  frequency  shift  Awjat,  linear 

frequency  wQ  and  growth  rate  y  versus  k/k^  for  Vg/v£.  *0.3.  We  note  that 
the  saturation  amplitude  ls4C  is  much  !•**  then  unity  only  if  k«  k^. 

However,  from  Eq.  (57)  we  observe  that  Aw.  .  approaches  -w  for  k«  k  . 

SIC  O  m 

This  violates  the  assumptions  of  our  perturbation  theory.  Nevertheless, 
it  is  true  that  |Wwl  is  largest  for  long  wavelength  modes  k«  k  , 
and  that  Eqs.  (50)  and  (55)  are  valid  only  for  je^/T.  | ,  jAw/wQj«1.  Also 
from  Eqs.  (37),  (39)  and  (40),  we  notice  that  wQ  becomes  smaller  than  the 


0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0  4.5 


K/K 


m 


Saturated  lower-hybrid  drift  mode  amplitude  le^/T;lsat  due  t0  non" 

linear  frequency  shift  mechanism,  normalized  frequency  tWu grow 

rate  y/u*a .  ,  and  nonlinear  frequency  shift  l iw  /id,.  as  functions  of 

k/k  for  v_/v  .  *0.3.  is  the  lower  hybrid  frequency,  and  k  is 

m  E  1 1  in  m 

the  wave  number  of  the  most  unstable  mode.  Note  that  jeb^/T.j  «1 
only  for  k/k  «1 ,  and  that  -Aw  _«  w_  only  for  k/k  >  1 . 
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Ion  cyclotron  frequency  for  k<0.4k  when  v_/v  .«0.3.  Therefore,  our 

-  m  E  1 1 

assumption  of  unmagnetized  ions  breaks  down,  and  these  long  wavelength 

lower-hybrid  drift  modes  convert  into  ion  cyclotron  drift  modes.  Many 

4-8 

studies  have  shown  that  a  nonlinear  ion  orbit  perturbation  induces 

a  frequency  shift  which  can  stabilize  the  ion  cyclotron  drift  instability. 

For  modes  with  k>k  ,  which  includes  the  most  unstable  mode,  k  ■  k  , 

-  m  m 

|e$^/T. |  ■  0(1),  which  is  too  large  and  thus  also  invalidates  our  per¬ 
turbation  theory.  For  these  short  wavelength  modes,  we  expect  that  other 
nonlinear  effects  will  be  the  dominant  saturation  mechanisms,  for  example, 
trapping  or  quasi  linear  diffusion. 

The  nonlinear  frequency  shift  does  not  appear  to  be  an  efficient 
mechanism  for  saturating  the  lower-hybrid  drift  instability.  However,  it 
is  important  to  compare  some  of  the  more  promising  saturation  mechanisms 
and  ascertain  whether  the  nonlinear  frequency  shift  effect  should  have 
been  Included  in  their  descriptions. 

Let  us  restrict  our  discussion  to  the  low  drift  regime  v_  «v  .  and 

E  ti 

T#«Tj.  By  using  Eqs.  (39),  (43)  and  (50),  the  field  energy,  £-<E^/8ir>, 
for  which  Aa)«-0.1w  (nominal  value)  is 


:  ■  i)(%  •  (^/ 


(8  et  +  5  v  v 

jit  x  m  m  'I  E  \ 
1  /  k2  \  Vti 


where  Aw  stands  for  the  nonlinear  frequency  shift  and  w»w  /(1+a  )  is  the 

o 

real  part  of  mode  frequency  including  the  linear  correction  term.  Note  that 
the  right  side  of  Eq.  (58)  is  a  function  of  k/k  and  v_/v.,,  and  not  a 

ffl  t  Cl 

function  of  T  /Tj . 
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An  oft  cited  saturation  mechanism  is  based  on  ion  trapping.  This 

1 2 

is  also  a  single-wave  effect.  For  v£>>vtj»  w'nske  et  al  estimated  that 
ion  trapping  requires 


where  C  is  order  of  0.1  to  fit  their  simulation  data.  For  the  low  drift 
regime,  i.e.,  Vg  < v  . ,  Chen  et  al1^  showed  that  above  the  threshold  given 
in  Eq.  (59),  saturation  occurs  when 


(  1  +  H)2  /u)2  ) 

pe  ce 


nT. 


i  / trap 


(60) 


This  estimate  was  calculated  on  the  basis  of  energy  conservation;  as  the 

ions  trap,  the  velocity  distribution  flattens  in  the  neighborhood  of  the 

wave  phase  velocity  ww/k  liberating  kinetic  energy  that  is  then  converted 
9  1 3 

into  wave  energy.  ’  Simulations  described  in  Ref.  13  demonstrated  that 
the  lower-hybrid  drift  instability  was  stabilized  by  ion  trapping  at  ampli¬ 
tude  consistent  with  Eq.  (60)  when  v^  was  kept  constant  in  time. 

Another  possible  saturation  mechanism  is  stabilization  via  current 

Q 

relaxation,  which  gives 
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This  mechanism  can  apply  to  single  wave  cases  or  turbulent  conditions. 


Finally,  saturation  due  to  electron  resonance  broadening 


10,11 ,14 


yields 


(  1  +  5 

pe  ce 


(62) 


for  plasmas  with  nonzero  .  This  equation  is  given  in  Ref.  11  without 

derivation.  This  mechanism  requires  the  presence  of  turbulence. 

Figure  3  is  a  plot  of  (fi/nT.)^  ,  ( £s/nTi )  trap  (£s/"Vcr 

versus  k/k  for  M/m  ■  3672  ,  v_/v  .  -0.3  ,  and  T  /T.-0.25  .  m  is 
m  t  1 1  e  i 

the  electron  mass.  It  is  seen,  except  for  k  «k  ,  that  (O/nT.)  is 

m  w  i  Aui 

much  larger  than  (P  /nT.)  and  (P  /nT.)  ;  and  thus  the  frequency 

3  '-'s  i  trap  &s  i  cr 

shift  is  important  only  for  k<0.1k  at  which  the  lower-hybrid  drift  mode  is 

m 

converted  into  a  drift  cyclotron  mode. 

In  Fig.  4  we  present  (£/"Vau,‘  (8s/"Vtrap  ,  (^/nT.)^ 

and  (Ss/nT[)rb  versus  vE/vt.  for  M/m » 3672  ,  T  /T.  -0.25  .  For  the 

parameters  chosen,  it  can  be  seen  that  for  v_/v  .  <0.2  the  lowest  satura- 

E  1 1 

tion  amplitudes  are  achieved  by  ion  trapping.  Again  it  shows  that  for  finite 

k/k  the  nonlinear  frequency  shift  is  insignificant  at  the  saturation  ampli- 
m 

tudes  suggested  by  the  other  proposed  saturation  mechanisms.  It  is  also 
noted  that  the  saturation  level  due  to  current  relaxation  is  only  slightly 
higher  than  that  of  electron  resonance  broadening  for  these  parameters; 
for  Te/T.  -0.25,  (£5/nT.)cri1.38(^/nT.)pb. 

In  conclusion,  we  have  shown  that  a  nonlinear  frequency  shift  has 
a  weakly  stabilizing  influence  on  the  lower-hybrid  drift  instability.  Ion 
trapping,  current  relaxation,  and  electron  resonance  broadening  are  more 
likely  saturation  mechanisms;  which  of  these  actually  accounts  for  satur¬ 
ation  depends  on  the  plasma  parameters. 
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,  e  p  /nT.  -U2/8^T.)  as 
.vbrid  drift  mode  amplitudes  b  *  ,q.25.  T*< 

ted  iower-hybn<J  .-0.3,  and  V  ,  .  a_d 

«  k/k  for  M/m- 3672,  V  tl  1aJtation  U-r.)  and 
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for  k/V°'  ,s  * 


Saturated  lower-hybrid  drift  mode  amplitudes  £■, /nT.  as  functions  of 
Vg/vt.  for  M/m  *3672  and  T^/T.  *0.25.  The  amplitudes  corresponding 
to  three  saturation  mechanisms  [ion  traping  (trap),  current  relaxa¬ 
tion  (c.r.),  and  electron  resonance  broadening  (r.b.)]  are  compared 
to  that  for  which  Au«-0.1ui  for  various  values  of  k.  Only  for 

k<0.1k  and  moderate  values  of  v_/v  .  are  the  frequency  shifts  signi 
m  t  1 1 

ficant  at  the  saturated  amplitudes  given. 
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